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Tj" Recently the HADES collaboration has published the invariant mass spectrum of 

o ■ 

' e+e pairs, dN/dMg+g-, produced in C+C collisions at 2 AGeV. Using electromag- 

<N : 

■ netic probes, one hopes to get in this experiment information on hadron properties 



at high density and temperature. Simulations show that firm conclusions on pos- 
sible in-medium modifications of meson properties will only be possible when the 
elementary meson production cross sections, especially in the pn channel, as well as 
production cross sections of baryonic resonances are better known. Presently one can 
I conclude that a) simulations overpredict by far the cross section at Mg+g- ~ M^^ 

• if free production cross sections are used and that b) the upper limit of the rj decay 

into e"*"e~ is smaller than the present upper limit of the Particle Data Group. This 
is the result of simulations using the Isospin Quantum Molecular Dynamics (IQMD) 
approach. 



I. INTRODUCTION 

Theory predicts since long that the properties of hadrons change if they are surrounded 
by matter. For baryons this change has been verified in 7A reactions where the total photon 



absorption cross section 



Jj shows a nontrivial dependence on the mass of the target nucleus. 



This nontrivial dependence has been interpreted as a change of the properties of the nuclear 



resonances in matter It is, however, difficult to asses whether the observed in- medium 
modifications have to be attributed to a change of the resonance properties or to a change of 
those of their decay products. Coupled channel calculations provide a mean to answer this 
question but presently neither the data are sufficiently precise nor the theoretical ingredients 
can be sufficiently well determined in order to allow for firm conclusions even if for some 
hadrons like the p and K 4] mesons a lot of progress has been made recently. 

The strategy is different for both cases. The study of the strange mesons takes advantage 
of the fact that they have to be produced in a heavy ion reaction, that each strange hadron is 
accompanied by an anti-strange one and that the production cross sections are phase space 
dominated. Systematic studies of the excitation function and of the system size dependence 
of the yields as well as of the modification of the measured K meson spectra as compared to 
that measured in pp collisions allow for conclusions on the interaction of the K's with the 
environment jsf. 

The p meson can decay into a dilepton pair which - being an electromagnetic probe 
- does not interact anymore with the nuclear environment. Therefore this dilepton pair 
carries direct information on the particle at the time point of its decay in the medium. 
The problem is that many resonances and mesons contribute to the dilepton yield and it is 
all but easy to determine which particle is at the origin of the dilepton pair. In order to 
compare data with theory, one has to identify all dilepton sources and their contribution to 
the dilepton spectra. This superposition of the different sources is called cocktail plot. If it 
deviates from experiment at least one of the sources is not correctly described and one may 
start to test how this source is modified by the hadronic environment. 

It was the DLS collaboration which first presented dilepton invariant mass spectra in 
heavy ion collisions at beam energies of around 1 AGeV [Si]- The systematic errors of these 
exploratory experiments have been, however, too large to allow for a detailed conclusion 
on the behavior of hadrons in matter. Later, at higher (SPS) energies, the CERES/NA45 
collaboration [3] presented spectra, which were not in agreement with the standard cocktail 
plots. Two theoretical models have been advanced to explain this difference. Rapp et al. ^ 
calculated the in medium modification of the spectral function of the p in hadronic matter. 
With this in-medium change of the spectral function the theoretical and experimental yields 
agree. As Eletsky et al. |8| explained p - meson and p - baryon interactions compensate 
each other as far as the shift of the pole mass is concerned but collisions broaden the 
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width considerably. Gallmeister et al. [oj] showed on the other side that the discrepancy 
disappears as well if one adds to the spectrum the emission of the dileptons from a thermal 
qq (or hadron-hadron) annihilation using lowest order QCD calculations. 

Most recently the NA60 collaboration measured very precisely the invariant mass spec- 
trum of dileptons in the p mass region [l^ but it is still debated whether the discrepancy 
between cocktail plot and data is due to a modification of hadronic properties or due to 



annihilation processes. Additional information may be obtained from the pt spectra [11 1 
because each emission source shows a specific transverse momentum pattern. However, con- 
sensus about the relative importance of the different possible production mechanism has not 
been obtained yet. 

To clarify this question it is necessary to study the dilepton production at lower energies 
where quarks remain bound in hadrons and hadron hadron annihilations are rare. Then the 
process proposed by Gallmeister is absent and thermal production does not play a decisive 
role. In addition one has to investigate small systems where direct collisions dominate over 
the production in the participant heat bath. 

Recently the HADES collaboration has published the dilepton invariant mass spectrum 



for the reaction C + C at 2 AGeV 



12|. This system is small and at this energy the formation 



of a quark phase is beyond reach as the analysis of many other observables has shown. It may 
therefore serve to solve the question of how the p meson changes in a hadronic environment 
provided that it can be proven that all the other ingredients of the cocktail plot are well 
under control. 

It is the purpose of this article to investigate in detail the dilepton invariant mass spectra 
using one of the presently available programs which simulate heavy ion reactions on an event 
by event basis, the Isospin Quantum Molecular dynamics (IQMD) approach. The main 
objective is to find out whether the present dilepton data are sufficiently precise to allow 
for conclusions on the theoretically predicted change of the particle properties in a nuclear 
environment or to identify the obstacles on the way to achieve this goal. We concentrate in 
this exploratory study on the most significant modifications: mass shifts and changes of the 
decay width. 

Before we present the results of our simulations we start out with a short presentation 
of the model and a discussion of our present theoretical and experimental knowledge on 
all the elementary processes which contribute to the dilepton spectra and of how they are 
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implemented in our simulation program. 



II. THE IQMD MODEL 

The semi-classical IQMD program 



13| simulates heavy ion reactions on a event by event 
basis and is one of the standard analyzing tools for heavy ion reactions at and below 2 
AGeV. In this program hadrons interact by potentials and by collisions. The former ones 
are Briickner G-matrix parameterizations for the baryons or parametrized meson-baryon 
potentials. Thus nuclei are bound objects with a binding energy following the Weizsacker 



mass formula. If two hadrons come closer than r = \J(7tot/'^ they collide. If several exit 
channels are available a random number determines which one is realized. The relative 
weight is given by the relative cross section. The momenta and the mass (if the particles 
have a finite width) of the hadrons in the final state are randomly determined. Their 
distribution follows either experimental measurements or phase space, if experimental results 
are not available. In the standard version [is!] of the program, nucleons as well as baryonic 
resonances, pions and kaons are the particles which are propagated. 

For the investigation presented here we have added production cross sections of all parti- 
cles which may contribute to the invariant mass spectrum of dileptons: np bremsstrahlung, 
7] Dalitz and direct decay, uo (Dalitz and direct) and p decay, A Dalitz decay and 7r° Dalitz 
decay. Because we concentrate on a very light system, where the probability that mesons 
have secondary interactions is small, it has not been necessary to add the (largely unknown) 
meson absorption or rescattering cross sections or to use off-shell transport approaches. 
When these particles are produced we use the branching ratios of the Particle Data Group 
to determine their contribution to the dilepton spectrum. 



III. ELEMENTARY DILEPTON CROSS SECTIONS 



A. vr production and decay 

1. -K^ decay into dileptons 

At low invariant mass the overwhelming number of dileptons comes from the decay of 7r° 
mesons which can decay into dileptons via 7r° — > e"^e~7. The shape of the mass distribution 
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of a dilepton in a vr" Dalitz decay is given by [l5|: 



m^ro is the mass of the ttq, me- the electron mass and M that of the dilepton pair. We take 
the branching ratio BR(7r° e^e~7) as 0.01198. 



B. t] production and decay 

In the energy regime which is of interest here, the rj production in pp collisions has been 
wen studied tKe TAPS Q a.d tKe DISTO co„abo.at,o. M. TK. can be seen in fig. [1] 
which shows on top the distribution of the rj excess energies in the nucleon nucleon collisions 
for the reaction C+C at 2 AGeV. The excess energy Xrj is defined as 

Xn = ^/s- 2Mn - M^. (2) 

We see that excess energies below 0.6 GeV are most relevant for this reaction. In the bottom 
)art of fig. [1] we display the world data points for rj production in elementary NN collisions 
IgI . ItI . 18|. Whereas the cross section a{pp pprj) is known over the whole excess energy 
interval which is relevant for our investigation the a{pn — > pnr]) cross section is known only 
up to an excess energy of a;^ = 0.12 GeV. Thus we have to extrapolate this cross section 
into the relevant excess energy domain. This extrapolation leaves a lot of freedom even 
if the rj meson production cross section has been measured in heavy ion reactions by the 
TAPS collaboration. The reason is that in heavy ion reactions a multitude of processes may 
modify the elementary cross section at the same nominal energy. These processes and the 
consequences will be discussed later. We parametrize the a{pn — > pnr]) and a{pp —>■ pprj) 
cross section by a fit using the form 

a{x.r,) = axj (3) 

with a = 1213.8, a = 162.1, a = 99.6 fib and b=1.50, b=-0.08, b=-1.24 for excess energies 
of Xrj < 283 MeV, 283 MeV < < 651 MeV , > 651 MeV for pp collisions and 
a = 25623, a = 324.3, a = 199 j^b and b=2.03, b=-0.08, b=-1.24 for excess energies of 
Xri < 200 MeV, 200 MeV < Xr^ < 651 MeV , > 651 MeV for np colhsions assuming that 
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FIG. 1: (Color online) Top: excess energy, Xr^, distribution of NN collisions in the reaction 2 AGeV 
C+C. Bottom: production cross section of the rj meson. The solid curves are fits (eq. [3]) to the 
data Q, Q 
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at large excess energies the np cross section is twice the pp cross section. These fits are also 
displayed in fig. [TJ We parametrize the shape of the mass distribution of the rj by 



o2 



with = 0.547 GeV and 1^=1.18 keV. 



1. Contribution of the N*(1535) 



The very detailed experimental investigation of the r] production in pp collisions at excess 
energies of 324, 412, and 554 MeV (corresponding to beam energies of i?beam=2.15, 2.5 
and 2.85 GeV) by the DISTO collaboration 17|] allows to identify the different production 
channels by analyzing the pr/ invariant mass spectrum. It turned out, as predicted by theory 
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22l | , that there are essentially two channels, a direct production channel and a production 



via the A^*(1535) resonance. The direct contribution follows the three body phase space for 
the pp — > ppr] reaction. The experimental mass distribution of the A^*(1535) resonance 
created in the reaction p p — A^*(1535) p can be described by a Breit-Wigner distribution 
of the form 17l |: 

<M) = , ,^.2^2^,.,. ^ (5) 

with 

x{M, Mr) = h^^^^^ji^ + K^^^ri^, (6) 



(M2 - M2)2 + MlT\x^{M, Mr) 



\^{Mr) ' """q^MRj 

where 6^ is the branching ratio of the decay N*(1535) —>■ N rj (which we assume to be 55%), 
b-^ is the branching ratio of the decay N*(1535) — N vr (which counts for 45%). and 
are the momenta of vr and rj in the frame of the resonance and are given by: 



and 



g^(M^.) 



\ 



Ml. 



- M"^ + M'^ 

2Mjv* ^ " 



Ml.-MS + M^^ 



2M, 



Ml 



N* 



(7) 



(8) 



We note in passing that in reference 17| the square on the x in eq. [5] has been forgotten. In 
: ig. [2]we display for the three energies which have been measured by the DISTO collaboration 
171] the total experimental and theoretical prj invariant mass distribution as well as the 
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different contributions to the theoretical curve. The experimental data are best reproduced 
for Mr = 1.530 GeV and F/j = 150 MeV. As expected, the iV*(1535) resonance enhances 
the low invariant mass part as compared to phase space. How the resonance production 





FIG. 2: (Color online) Simulated invariant mass spectrum of the outgoing proton and the rj meson 
in the p p — > p p reaction for three beam energies (2.15 GeV, 2.5 GeV and 2.85 GeV). The curves 
represent the sum of the contributions from the two productions channels of the rj, the direct 
production and that via the N*(1535) resonance. The data [17^ have no absolute normalization. 
We normalize them here to our result at M7v*(i535) = 1500 MeV/c^ 

modifies the spectra of proton and rj as compared to the production according to the three 
body phase space is shown in fig. [31 On the left hand side we display the center of mass 
momentum of the t], on the right hand side the proton momentum in the pp rest frame. 
Choosing these variables allows for a comparison with the experimental results. We see 
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clearly the consequence of the t] resonance production and therefore it will be difficult to 
separate the modification of the 1] in the medium from that of the N*(1535) resonance. Both 
will show up as a modification of the dilepton spectra. 



2. T] decay into dileptons 



With a branching ratio of 6 ■ 10 ^ iJ] the rj decays into e+e 7. The shape of the invariant 
mass distribution of the dilepton pair is given by [15| 



dN 

dM 



-(1 + 2^ 



M2 

)(!- — 
mi 



1-4 



M2 



(9) 



is the mass of the rj, rrie- the electron mass and M that of the dilepton pair. It has been 
shown that this expression has to be multiplied with an electromagnetic form factor. With 



dM 



dN 
dM 



where 



F{M^ 



-1 _ Ml' 

A'i 



(10) 



(11) 



with = (0.72 ± 0.09) GeV one finds good agreement with data [20|. In addition to the 



three body decay there may also be a two body one into a dilepton pair. The Particle Data 
Group ij] quotes as an upper limit a branching ratio of 7.7 ■ 10"^. We include this value in 
our standard calculation (standard will be explained later). 



C. uj production and decay 



The LJ production in pp collisions for excess ener gies below 440 MeV has been studied at 
COSY 23!, at SATURNE ^ and by the DISTO 251 collaboration. The cross section as 
well as our fit of the form ax^j'^ where x^j is the excess energy in MeV, a = (192.204 ± 8.622) 
/i6 and b = 1.12182 ± 0.1077 is shown in fig. HI We include in our simulation as well the 
endothermic {^/sq ^m^-m.n-= 643 MeV) reaction vr+N^ cj+N. Because vr's have usually only 
a small energy this reaction is less important than the baryonic channel. The experimental 



data have been parametrized 



26| by 



cr-KN-.ujNi.rnh) 



1.38(yi 



1.6 



0.0011 + (x/i-^) 



1.7 



(12) 
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FIG. 3: (Color online) Differential cross section in IQMD as a function of the center of mass 
momentum of the ij meson (left) and as a function of the proton momentum in the pp rest system 
(right) for pp — > ppT] collisions at different beam energies, Ebeam = 2.15 GeV (up), Eheam = 2.5 
GeV (middle) and Ei^^am = 2.85 GeV (down). Solid lines represent r/ production including the 
contribution of the A^*(1535) resonance and dashed curves represent the direct production via an 
uniform three body phase space distribution. The experimental data are form ref. [l3] 
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with a in mb, ^/s and ^/sq in GeV. ^/sq = tun + rn^ (1.721 GeV for uj in vacuum) is 
the threshold energy. It has been suggested that the production of the oj passes by the 



excitation of baryon resonances 



27 



where the A^" 



29|, 



535) plays a prominent role having 



3G, 



3l| . It produces uj mesons with 



a substantial branching ratio into the Nuo channel 
masses well below 783 MeV. If this were the case the strong cjN coupling would lead to a 
strong off-shell contribution to da/dM (M being the invariant mass of the dilepton pair) at 
invariant masses well below the free u mass peak. This off-shell uj production would even 
dominate the dilepton spectra up to excess energies of several hundred MeV. Only very 
recently calculations of the spectral function have been advanced which exploit the available 



30 



32|. 



7A^ and vrA^ data in a coupled channel analysis 

A J A — > u) experiment was recently performed by the CBELSA/TAPS collabora- 
tion. They observed that the pole mass decreases with increasing density of the envi- 
ronment. 33|. For momenta less than 500 MeV/c^, they observed an u pole mass of 
M=[722li(stat)lf (syst)] MeV/c^ for an average density of 0.6 po- Unfortunately, no sig- 
nificant measurement of the width was obtained due to the dominance of the experimental 
resolution. Using this data and the Brown-Rho scaling formula: 



* (1 - a 



Po' 



(13) 



we find a=0.13. Fig. [5] shows the density distribution at the uj production points for a C+C 
collision at 2AGeV. The average density of < p >= 1.394 po is twice as large as for the TAPS 
experiment. Applying eq. [2] we obtain a wide distribution around the average pole mass of 
M= 641 MeV. 

In our simulation we have the option to use this in-medium mass modification. Because 
there are no conclusive results on the width we kept the free value of 8 MeV. The shape of 
the invariant mass distribution of dileptons from the uj decay is given by the Breit-Wigner 
distribution : 



dN 
dM 



:i + 2- 



Af2 



M2 



{ml - M2)2 + [m, 
with T^='& MeV and m^^ as defined in eq. fT3l 



(M2/4-m2_)3/2^^2 



(14) 



M (m2/4-m2_)3/; 



Another uncertainty is the production of the uj in pn reactions. In meson exchange 
models the relative strength of the production in pp and pn reactions depends strongly on 
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FIG. 4: (Color online) Production cross section of the lo in pp co^ 
440 MeV and our fit of the form a = axf*,. The data are from ref 



lision up to an excess energy of 



23 



24 




1.2 1.6 2.0 



FIG. 5: (Color online) Distribution of the density at the omega production points in units of the 
normal nuclear matter density po for the reaction C+C at 2 AGeV. 

the quantum number of the exchanged mesons. Neglecting possible differences due to initial 
and final state interactions, we expect a{pn — * pnuj)/a{pp —>■ ppu) = 5, if only isovector 
mesons (vr, p) are exchanged 3J]. The two data points for the reaction np du point toward 
an enhancement of the pn cross section as compared to the pp cross section ^]. The error 
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bars are, however, too large in order to quantify this enhancement. In our simulations we 
assume a{pn —>■ pnuj) = b* a{pp ppuj) with different values of b. 

The u contributes to the dilepton spectrum in two different ways. Either it decays directly 
into a dilepton pair whose invariant mass equals that of the u meson or the dilepton pair is 
accompanied by a ttq meson. For the latter channel the shape of the dilepton invariant mass 



distribution has been parametrized by KroU 15(] 



^ = ^(1 + + ? - 4 l - 4^ (15) 

with M being the invariant dilepton mass. This Dalitz type decay has to be corrected by an 



electromagnetic form factor 20 1 



^dM' ' ' dM 

with 



'tot = Fm * — (16) 



and a = 0.6519 GeV , b = 0.04198 GeV in order to be in agreement with data. The 
branching ratios into the two channels are given by 5.9 ■ 10^^(7.14 ■ 10^^) for the e"'"e~7r 
(e"'"e~) channel jl^. Both, the unknown pn cross section as well as the little known off- 
shell contribution at small excess energies make it difficult to predict the io contribution at 
invariant dilepton masses between 0.6 and 0.8 GeV. 

D. p production and decay 

In our simulation the p meson can be produced in three channels: A^A^ — > NNp, ttN — »• 
pN and tt+tt^ p. 

The few experimental data points of the total cross section in the A^A^ —>■ NNp channel 
have been fitted by jssl 



0.24(v/i-^) 

with y/so = 2.646 GeV being the threshold of the reaction. In view of the strong coupling 
of the p to nuclear resonances this course-grained parametrization has most probably large 
systematic errors and presents a lower limit to the p production. Other models like URQMD 
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use a parametrization of the resonance production which yield higher p yields. For the tt+N 



N+p data [s^ we use the parametrization of [26 1 

1.5(v/i 



(^■KN^pN{rnh) 



2.2 



i3.5 



(19) 



0.0018 + (v^- V^oJ 
with ^ = 1.708 GeV. 

Having a large width and therefore a short life time, the p meson is an ideal particle to 
probe whether the nuclear environment changes mesonic properties. If produced in hadronic 
matter the majority of them decay in matter and therefore the dileptons carry direct in- 
formation on the in-medium properties. Theory predicts that these properties are different 



De now consensus that the width of the 
3, 
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37l |. the question of how the pole 



from that of the free p. Whereas there seems to 
p increases if brought into a nuclear environment 
mass changes is still debated. Based on QCD sum rule calculations, Hatsuda and Lee [38 1 
predicted a lowering of the p mass in a nuclear environment, a suggestion which has later 



been confirmed by Brown and Rho |39|, 



More recent and more sophisticated calcula- 



371 . Experimentally the 



tions leave, on the contrary, the p mass almost unchanged 
situation is also far from being clear. In pA collisions [41] at 12 GeV a decrease of the mass 
(m(p)/m(0) = 1 — 0.09p/po - about half of the value predicted by theory) and no increase 
of the width has been reported. The dilepton data in In+In collisions at 158 AGeV [l^ 
are best described using the free p pole mass but a considerable broadening of the mass 
distribution. In contradiction to the earlier theoretical expectations this broadening is al- 
most symmetric around the pole mass but recently it has been pointed out 4^ that the 
^-functional approach may explain this symmetry. 

Whether these experimental differences are exclusively due to the different environments 
(cold nuclear matter in pA reactions, an expanding meson dominated fireball after a possible 
phase transition from a quark gluon plasma in AA collisions) has not been fully explored 
yet. It is very difficult to exploit this experimental information for heavy ion reactions at 
2 AGeV where theory predicts that most of the p mesons are decay products from nuclear 
resonances, especially of the A^*(1520) resonance which has a branching ratio of 15 — 25% 
into the pN channel. For the present status of the theoretical spectral function calculations 
for the p meson we refer to 32, l43[ |. 

As for the uj meson the inconclusive situation of theory and experiment suggest to employ 
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for this exploratory study the free pole mass distribution of the p: 



dN 
dM 



(20) 



with Trip = 0.775 GeV/c^ 



m 



P2 



0.761 GeV/c2, Tp = 0.118 GeV/c^ ^ and the 



parametrized free cross sections. For the branching ratio of the p into dileptons we use 
4.5 ■ 10-^ 



E. pn - bremsstrahlung 

In each np collisions real and virtual photons can be produced. The invariant mass 
distribution of the e~^e~ pairs, the decay product of the virtual photon, is given by: 

dP{s,M) _ laM s-(mp + m„)2 .Qmax + qo max Qmax \ /m \ 

dM ^ 3^M M ^ ^ 

with 

_ g + - (mp + TJlnf f . 

QOmax — n /~ \ ) 

Za/ S 



'qlmax - M'. (23) 

^/s is the np center of mass energy, Ccm is the energy of the incoming proton in the np center 
of mass system, a is the electromagnetic coupling constant, rrip and m„ the masses of proton 
and neutron, qomax the maximal dilepton energy and Qmax the maximal dilepton momen- 
tum. The bremsstrahlung from pp collisions is of quadrupole type and can be neglected as 
compared to the dipole pn bremsstrahlung. 

F. A Dalitz decay 

It is not experimentally verified yet whether the Dalitz decay into e~^e~ of the A resonance 
exists but since it decays into a photon it should also decay into a dilepton. The width of 



the Dalitz-decay to dileptons of invariant mass M is determined by QED 45| : 

dT _ a ro(M^) 

Tap ~ 3 TT M2 

where 



(24) 



T,iM') = ^'^'^f5'''"T'"^^ ^^[2 mM') + M,iM')] (25) 
lo TT 
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is the total decay rate into a virtual photon with mass M and 

A(a;, z) = + + — 2{xy + xz + yz). (26) 

Mt and Mi depend on the form of the interaction. For the A decay we take the A^A7 vertex 

n 

from Using this interaction we obtain the following matrix elements: 

2 

Ml = {ef gf^^M^ A{m^ - - go) 

2 

Mt = {ef gf^P^[qo\5 - 3(go + m^)) - M\m^ + + go)] (27) 

with 

f = -1 . mA + mjv 

^ m^((mjv + mA)2-M2)' ^ ^ 

go the energy of the dilepton pair in the A center of mass, e the electric charge and 
g = 2.72 is the coupling constant fitted to the photonic decay width ro(0) = 0.72 MeV 

IV. THE C+C REACTION AT 2 AND 1 AGEV 

For the simulation of the heavy ion reaction we use the IQMD prog ram which has been 
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described in section I. The details of this program can be found in 

The presented results are impact parameter averaged and have been corrected for the 
experimental mass resolution and acceptance with a program provided by the HADES col- 
laboration. We have neglected in our calculation the reabsorption cross section of the r) 

i I 

mesons which is of the order of 20 mb [48j] in our kinematic domain but of little importance 
for such a light system. We compare the results of the standard set up, where free masses 
and widths as well the most common extrapolations or theoretical predictions of unknown 
cross sections are used, with calculation in which it is assumed that the particle properties 
change in the medium or in which other cross section parametrization are applied. 

Fig. [6] shows the result of the standard simulation set up: a{np — > nprj) = 2cr{pp ppt]) , 
a{np — > npuj) = 5a{pp ppuj), = M° and the branching ratio BR^^e+e- = 7.7- 10~^. It 
is called standard because it uses standard literature values for the unknown physical input 
quantities. We see first of all that with the resolution of the HADES experiment the direct 
7] decay would yield a visible peak which is not present in the data. Therefore the upper 
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limit has to be lower than that quoted by the Particle Data Group [l^ . We see as well that 
the simulations overpredict the yield in the region of the free uj mass. This confirms the 
result of the simulations with other programs which have been published by the HADES 



collaboration 



121]. On the contrary, the simulations reproduce well the mass region in which 



the lepton pairs are coming dominantly from the rj decay and from pn bremsstrahlung. If 
the experimentally unknown a{np pnu) equals a{pp ppuj) the yield in the Ai^ mass 
region would be strongly reduced and comes closer to experimental data as can be seen 
in fig. [71 We obtain the same level of agreement with data if we take a{np pnuo) = 
5(j{pp — > ppu) but assume in addition that the mass of the u decreases in the medium 
according to eq. [131 as indicated by the CBELSA/TAPS results |33|]. This can be seen in 
fig. [HI The experimental error bars are large, however, do not show much structure and 
the deviations are only a factor of two. Higher statistics data would certainly improve this 
situation. The best agreement is obtained in simulations with a{np pnuo) = a{pp ppuj) 
and an in medium uo mass, as seen in fig. [9l 

Therefore, without further informations on (j{np — > npuj) heavy ion reactions will not 
reveal any robust information on in medium modifications of the uo meson. 

If we assume a{np —>■ nprj) = a{pp pprj) in the region where no data on the a{np — *■ 
npt]) cross section is available (only 5% of the rj are produced at an energy where experimental 
information on this cross section is available), we underpredict slightly the yield in this mass 
region as seen in fig. [lOl The experimental error bars are too large, however, in order to 
conclude more than that there are indications that if the mass of the rj does not change in 
the medium the a{np —>■ nprf) is larger than (j{pp — > pprf) at excess energies above 100 MeV. 
As for the uj meson possible in medium changes of the rj meson require a detailed study of 
its production in the pn channel. 

Fig im summarizes the study of the influence of the parametrization of unknown processes 
on the dilepton yield. If one compares the results of the different scenarios of table HI with 
the experimental results, we see that the standard parameterizations of these input quantities 
(A) yield not results which are in agreement with data at invariant dilepton masses around 
550 MeV and in between 750 MeV and 950 MeV. The former difference suggests that the 
partial width for the disintegration of the r] into a dilepton pair is much smaller than the 
upper limit quoted by the Particle Data Group [1^. The latter discrepancy contains the 
interesting physics as far as in medium particle properties are concerned. We see that even 
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FIG. 6: (Color online) The invariant mass spectrum of the HADES collaboration as compared 
with IQMD simulations for C+C at 2AGeV using a{np npr]) = 2a{pp pprf) , a{np — > npto) = 
5a{pp — > io), Mu = and the branching ratio {r] — > e+e") = 7.7 10~^ (model A). 




FIG. 7: (Color online) The invariant mass spectrum of the HADES collaboration as compared 
with IQMD simulations for C+C at 2AGeV using a{np npr]) = 2a{pp pprf) , a{np — > npuj) = 
a{pp ppto) = and the branching ratio (ry e^^e") = 7.7 10^^ (model B). 

a reduced uj production cross section in the np channel (B) does not render the calculation 
compatible with the data. Also the assumption that the mass of the lj changes in the medium 
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FIG. 8: (Color online) The invariant mass spectrum of the HADES collaboration as compared 
with IQMD simulations for C+C at 2AGeV using a{np npr]) = 2a{pp pprf) , a{np — > npto) = 
5a{pp ppuj) = M°(l — 0.13p/po) and the branching ratio (77 e+e~) = 7.7 10~^ (model C). 
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FIG. 9: (Color online) The invariant mass spectrum of the HADES collaboration as compared 
with IQMD simulations for C+C at 2AGeV using a{np npr]) = 2a{pp pprf) , a{np — > npu) = 
a{pp ppuj), = M^{1 — 0.13/9//9o) and the branching ratio (77 e'^'e") = 7.7 10~^ (model E). 

but that it is produced with the free cross section (C) overpredicts the experimental results 
because it shifts the surplus only to lower invariant masses. Only the combination of a 
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FIG. 10: (Color online) The invariant mass spectrum of the HADES collaboration as compared 
with IQMD simulations for C+C at 2AGeV using a{np nprf) = a{pp — > ppr]) , a{np — > npuj) = 
5a{pp — > ppuj) = M^(l — Q.12>p/ po) and the branching ratio (ry — > e+e^) = 7.7 10^^ (model 
D). 

lower in medium mass and a reduction of the standard assumption on the cross section in 
the pn channel (E) yield results which are compatible with the experimental error bars. 
The scenario (D) demonstrates that the data are not sufficiently precise to allow for robust 
conclusions on the np npr] channel. A variation of a factor of two gives results which are 
both compatible with experiment. 

Thus the C+C data at 2 AGeV show interesting new physics which is not compatible with 
the input of state of the art transport codes. Unfortunately without further information on 
the elementary cross sections with a neutron in the entrance channel it will not be possible 
to identify the origin of this discrepancy because a modification of the mass of the mesons 
in the medium yields the same effect as a change of the (experimentally unknown) cross 
section in the np channel. 

Lowering the energy to 1 AGeV the importance of the different channels changes and a 
comparison between the 2 and 1 AGeV data will elucidate part of the physics. Because the 
experimental data are divided by the number of vr^ the spectra for the pions change only 
little due to the acceptance corrections. The same is true for the A Dalitz decay. The yield 
of e~^e~ pairs from t] Dalitz decay and bremsstrahlung are lower, on the contrary, and the lo 
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FIG. 11: (Color online) Invariant mass dependence of the ratio theory/experiment for the C+C 
reaction at 2 AGeV for the different parametrization of unknown physical input quantities (see 
table III for details) 

production is practically absent due to the lack of energy (even if one takes into account that 
the Fermi momentum may create a larger ^/s value than in NN collisions at the same beam 
energy). Fig. [12] displays our filtered and acceptance corrected results. In the intermediate 
mass region the A Dalitz decay and bremsstrahlung have gained importance and are of the 
same order of magnitude. Dilepton pairs from t] Dalitz decay are less frequent and are not 
dominant anymore in the intermediate mass region. At this energy about 35% of the r) 
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come from a ^/s region where the np production cross section is known. So the uncertainty 
in this channel is reduced but still present. In the standard set up of the simulations the 
dilepton invariant mass spectrum at intermediate masses has always a strong component of 
the bremsstrahlung which gives about 50% of the yield. At invariant masses of around 200 
MeV the A Dalitz decay contributes the other 50% - if it exists. The data at 1 AGeV should 
therefore allow do define an upper limit of the A Dalitz decay. At higher invariant masses it is 
the T] decay which contributes the other 50%. If we assume that a{np —>■ nprj) = a{pp ppv) 
the 7] yield becomes that low that its influence on the spectrum is hardly visible. If we lower 




FIG. 12: (Color online) The invariant mass spectrum of IQMD simulations for C+C at lAGeV 
using a{np nprj) = 2a{pp ppif) , (T{np — > npu) = 5a{pp PP^), = and the branching 
ratio (r/ e+e~) = 7.7 10"^. 

the in medium u mass we see a larger u production cross section but it remains a small 
contribution to the total yield, as seen in fig. [131 

It is interesting to see in detail the differences between elementary collisions at ^/s = 
2.697 GeV and heavy ion collisions at the same nominal energy which show the large y/s 
distribution of the NN collisions displayed in fig. [TH There we see two peaks. The high 
energy peak is due to collisions between projectile and target nucleons, whereas the low 
energy peak is due to collisions among either projectile or target nucleons. The latter 
collisions contribute only to the bremsstrahlung and to the vr'' part of the dilepton spectrum. 
Due to rescattering the maximum of the distribution of the primary collisions is shifted 
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FIG. 13: (Color online) The invariant mass spectrum of IQMD simulations for C+C at lAGeV 
using a{np — > np-q) = 2a{pp ppr]) , a{np npto) = 5a{pp — > ppuj), = M^{1 — 0.13p/po) and 
the branching ratio (t/ — > e~^e~) = 7.7 10~^. 

toward a lower ^/s value. The consequences of the broad ^/s distribution on the t] and lo 
production as compared to elementary collisions at the nominal energy are summarized in 
table [H The first line shows the average ^/s value of all collisions above threshold. For the rj 
this value is slightly below, for the u - due to the larger threshold - slightly above the value 
for elementary collisions. Also the average number of collisions in C+C reactions depends 
on the particle type as seen in the second line. For the rj production we find 4.65 collisions 
above threshold for the u production 2.32. For the standard scenario {m^ = m°, BR(r7 — »• 
e~^e~) = 7.7.10^^, cr(pn pnif) = 2a{pp — ppr]) and a{pn — > pnu) = 5a{pp ppu)), we 
display in the third and the fourth line the average production cross section in np and pp 
collisions in the heavy ion reaction as compared to the elementary reaction. For the 1] the 
average a{pp — > pprj) and a{pn —>■ pnt]) are lower in CC collisions than in elementary ones. 
This decrease has two origins: firstly the lower < \fs coii>threshoid > ^^"^ secondly the form of 
the ?7 production cross section which has a maximum at around y^=2.697 GeV and stays 
almost constant at higher energies. For the lo meson the situation is completely different. 
The elementary cross section increases with energy for all relevant energies and the average 
value in C+C is larger than that in elementary collisions. Therefore np as well as pp 
collisions in the heavy ion reaction produce more lo mesons than elementary collisions at the 
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same nominal energy. Consequently the enhancement factor of rj and uj mesons in heavy ion 
collisions is very different. 

For the 1 AGeV reaction (table [IT!) the situation is very different. The ^/s distribution of 
the collision in C+C is displayed in fig. [151 In elementary NN collisions at the same energy 
neither uj nor r] mesons can be produced ( threshold ^ = 2-659 GeV and A^.e.hoid „ = 
2.424 GeV). However, with the Fermi momentum, in C+C collisions subthreshold lo and r] 
production is possible. Due to the larger threshold uj production is suppressed with respect 
to the ?7 production. The production cross section at this energy tests the Fermi motion 
in the simulations which is not easy to model in semi-classical simulation codes. Therefore 
systematic errors reduce the predictive power for the meson production at this energy, but 
the analysis of the subthreshold kaon production shows that in between a factor of two the 
results are certainly trustworthy. 

In summary we have shown that the dilepton spectrum measured by the HADES collab- 
oration in the reaction 2 AGeV C+C at invariant masses above 600 MeV is not compatible 
with the standard scenario of simulation programs which uses free cross sections and free 
meson masses. Introducing a medium modification of the uj mass and lowering the unknown 
'pn —>■ pnuj cross section brings the calculation in agreement with data. The extrapolation 
from elementary cross section at the same nominal energy to heavy ion reactions is all but 
trivial. It depends on the threshold and on the energy dependence of the cross section. 
Before the elementary production cross sections in pn reactions are not determined and be- 
fore the cross sections for baryonic resonances are not better known heavy ion data do not 
provide the desired information on possible in medium modification of the meson properties. 
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